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Self-consistent random phase approximation (RPA) approaches in the relativistic framework are 
applied to calculate the isospin symmetry-breaking corrections Sc for the 0"*" 0"*" superallowed 
transitions. It is found that the corrections Sc are sensitive to the proper treatments of the Coulomb 
mean field, but not so much to specific effective interactions. With these corrections Sc, the nucleus- 
independent J^t values are obtained in combination with the experimental ft values in the most 
recent survey and the improved radiative corrections. It is found that the constancy of the Tt values 
is satisfied for all effective interactions employed. Furthermore, the element Vud and unitarity of the 
Cabibbo-Kobayashi-Maskawa matrix are discussed. 
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I. INTRODUCTION 

The Cabibbo-Kobayashi-Maskawa (CKM) matrix 
0] relates the quark eigenstates of the weak interaction 
with the quark mass eigenstates. The unitarity condition 
of the CKM matrix provides a rigorous test for the stan- 
dard model description of electroweak interactions. Its 
leading matrix element, Vud, only depends on the first 
generation quarks and so is the element that can be deter- 
mined most precisely. There are three traditional meth- 
ods to determine \Vud\ experimentally: nuclear 0"^ —f 
superallowed Fermi (3 decays [1], neutron decay (3], and 
pion P decay Q. Recently, experiments with nuclear 
mirror transitions provided another independent sensi- 
tive source for extracting the value of \Vud\ 0- 

Among these methods, the most precise determination 
of \Vud\ comes from the study of nuclear O"*" — > 0+ super- 
allowed Fermi /? decays Q . These pure Fermi transitions 
between nuclear isobaric analog states (IAS) allow for a 
direct measurement of the vector coupling constant Gy 
of semileptonic weak interactions by 

Together with the Fermi coupling constant Gf for purely 
leptonic decays, the up-down element of the CKM ma- 
trix can be determined, Vud = Gv/Gp- In Eq. ([T|), 
Kl{hcf = 27r3?iln2/(meC^)^ and is the transition- 
independent part of radiative corrections caused, for ex- 
ample, by the processes where the emitted electron may 
emit a bremsstrahlung photon that goes undetected in 
the experiment Q . The nucleus-independent Tt value 
is obtained by the corrections to the experimental /i 
values for radiative effects as well as isospin symme- 
try breaking by Coulomb and charge-dependent nuclear 
forces 0, 

.^i = /i(l + 5^)(l + <5NS-'5,), (2) 



where / and t represent the statistical rate function and 
partial half-life, respectively, and are obtained through 
measurements of the Q values, branching ratios, and half- 
lives for the superallowed decays. The correction terms 
(5^ and 5ns represent the transition-dependent radiative 
corrections Q . The correction term 8^ is the isospin 
symmetry-breaking correction, accounting for the isospin 
symmetry breaking in nuclei. 

The isospin is not an exact symmetry mainly due to 
the presence of the Coulomb forces in nuclei. The non- 
conservation of isospin symmetry induces a slight reduc- 
tion of the superallowed transition strength |A/i?p from 
its ideal value |Mop: 

l^/Fp = |(/|T±|^)p=.|Mop(l-5e), (3) 

where Afo = \/2 for T = 1 states with the exact isospin 
symmetry. 

Shell model calculations are generally used to deter- 
mine the isospin symmetry-breaking corrections be- Re- 
cently, by including the core orbitals, an improvement 
on such corrections has been achieved and a good agree- 
ment among the nucleus-independent Tt values for the 
13 well-measured cases has been obtained 0. 

Alternatively, the self-consistent random phase ap- 
proximation (RPA) based on microscopic mean field the- 
ories is another reliable approach for the superallowed 
transition strength Mp. Such calculations were per- 
formed for a few nuclei with the non-relativistic Skyrme 
Hartree-Fock approach in the 1990s [l3|. Since then 
no further investigation followed even though significant 
progress in self-consistent RPA in chaxge-exchange chan- 
nels have been made [ll|, ESj Elj H, [l3 ■ 

During the last decade, great efforts have been dedi- 
cated to developing the charge-exchange (Q)RPA within 
the relativistic framework. From the early model which 
only contains a rather small configuration space to 
the sophisticated model which includes Bogoliubov trans- 
formation and proton-neutron pairing [14| . these ap- 
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proaches are aimed at describing the spin-isospin reso- 
nances, f3 decay rates, neutrino-nucleus cross sections, 
etc., in a systematical, reliable, and predictive way. Re- 
cently, based on the success of the newly established 
density -dep endent relativistic Hartree-Fock (RHF) ap- 
proach [la 113, a fully self-consistent charge- exchange 
RPA has been established and the first applications have 
been performed for spin-isospin resonances like Gamow- 
Tcller and spin-dipolc resonances [15[- A very satis- 



factory agreement with the experimental data was ob- 
tained without any readjustment of the energy func- 
tional. Therefore, it is appropriate now to re-investigate 
the isospin corrections for superallowed Fermi /3 decay 
with these relativistic approaches. It is not the aim here 
to claim that a covariant framework is necessarily more 
appropriate for this problem than a non-rclativistic one 
such as Skyrme Hartrce-Fock plus RPA. The key point, 
which will be discussed in Sec. IIII A| is a full treatment 
of Coulomb and nuclear interactions in both their direct 
and exchange contributions. In this respect, satisfactory 
non-relativistic RPA studies of the Sc corrections are not 
available. 

In this paper, the self-consistent RPA approaches in 
the relativistic framework will be applied to calculate 
the isospin symmetry-breaking corrections Sc- With the 
corrections thus obtained, the nucleus-independent J-t 
values will be obtained in combination with the experi- 
mental ft values in the most recent survey (loj and the 
improved radiative corrections 0, 01- The element Vud 
and unitarity of the CKM matrix will then be discussed. 



II. SELF-CONSISTENT RELATIVISTIC RPA 

The basic ansatz of the relativistic Hartree (RH) , also 
known as relativistic mean field (RMF), and relativis- 
tic Hartree-Fock (RHF) theories is a Lagrangian density 
C, where nucleons are described as Dirac spinors that 
interact via the exchange of cr-, w-, p-, 7r-mesons and 
the photon [2l|, H^l ■ In order to give a satisfactory 
description of nuclear matter and finite nuclei, the non- 
linear self-coupling of mesons, e.g., in Refs. psi. \2^. [25j . 
or density-dependent meson-nucleon couplings, e.g., in 
Refs. [H, nil, are introduced. 

The effective Hamiltonian operator H can be obtained 
with the general Lcgendre transformation. Together with 
the trial ground state (Slater determinant) in the Hartree 
or Hartree-Fock approximation, the energy functional 
can be written as 



$0 Iff 



$0 



a ab 

-I^(a6|ni,2)|a6), 



(4) 



where the first term is the kinetic energy, the second 



and the last terms are the direct (Hartree) and exchange 
(Fock) energies, respectively. In the Hartree approxima- 
tion, the Fock term is neglected for simplicity. The two- 
body interaction V{1,2) includes the following meson- 
nucleon and photon-nuclcon interactions: 



K(l,2) 

Vp{l,2) 
^.(1,2) 



^a(1,2) 



= -[5<T7o]i[5<T7o]2-Da(l,2), 

[gu;707'"]l[guj707tJ.]2Du,{l,2), 

[5^707^7^]! ■ [gp7o7M^]2£'p(l,2), 



(5a) 
(5b) 
(5c) 



^.(1,2), 



A 

^707^1 -^3)]i[7o7m(1 -^3)]2i?A(l,2) 



(5d) 



with the finite-range Yukawa type propagator 

2 g-mi|ri-r2| 



A(l,2) 



An |ri-r2| 



(5e) 



(6) 



Furthermore, in order to cancel the contact interaction 
coming from the pion pseudovector coup ling , a zero-range 
pionic countcrtcrm should be included |15l. |21|: 



V^s{l,2) 
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(7) 

with g' = 1/3. Thus, g' is not an adjustable parameter. 

The RPA equations can be obtained by taking the sec- 
ond derivative of the energy functional E. In the charge- 
exchange channels, the RPA equations become 



'^ripn'p' 



p'n' 



^ pn 
*^ np 



(8) 



where p and p (n and n) denote unoccupied and occupied 
proton (neutron) states. These equations describe both 
the r+ and T_ channels. It should be emphasized that 
the unoccupied states include not only the states above 
the Fermi surface, but also the states in the Dirac sea. 
The RPA matrices A and B read 

^2,34 = (^1 - £;2)<5i2,34 + (14| Vph |32 - 23) , (9a) 
Si2,34 = - (13| 1/ph |42 - 24) , (9b) 

where the first term in the ket represents the direct con- 
tribution, and the second term represents the exchange 
contribution. In the RPA built on the Hartree mean field, 
the exchange contributions in Eqs. ([9]) arc accordingly 
neglected. 

In the self-consistent RPA calculations, the particle- 
hole residual interaction Vph should be derived from the 
same energy functional E as that used in the ground- 
state description. The explicit density dependence of the 
mcson-nuclcon couplings introduces, in principle, addi- 
tional rearrangement terms in the particle-hole residual 
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interaction Vph, and their contributions are essential for 
a quantitative description of excited states [13 ■ How- 
ever, since the rearrangement terms are due to the de- 
pendence on isoscalar ground-state densities, it is easy to 
see that they are absent in the charge- exchange channels. 
Therefore, in the description of superallowed Fermi /? de- 
cays, the particlc-holc residual interaction Vph is just the 
meson-nucleon interactions shown in Eqs. (|5ap - (|5dp and 
(O . The photon- nucleon interaction in Eq. ((5c)) does not 
contribute to the particle-hole residual interaction be- 
cause the configurations are of the neutron-proton type. 

The eigenvectors of the RPA equations ^ are sepa- 
rated into two groups, which respectively represent the 
excitations of the T_ and T+ channels with the following 
normalization conditions: 



Epn(C^/fi )J - EnpiK^l = +1- foi- channel, 



-1, for T-f channel. 



(10) 

Then, the excitation energies and X, Y amplitudes in 
the T_ channel read 



^ pn 1 np 



V, 



Jv 



(11) 



whereas the excitation energies and X, Y amplitudes in 
the r_(. channel arc 



np np 1 



1 u . 



(12) 



The 0^ — > 0"'' superallowed transition operators arc T_ 
or T-f . The transition probabilities between the ground- 
state and excited states read 



B 



Jv 



BSv- 



E 

p7l 



+ > (-)^ 



np 



(13a) 



^^{n\\T+\\p) 



(13b) 



Before ending this section, it is worthwhile to make 
the following remark about the self-consistency of the 
RH+RPA approach when it is applied to the 0"'' 0+ 
transitions. Within this approach, it is known that, in 
order to reproduce the excitation energies of Gamow- 
Tcller resonances, one has to adjust the ttNN particle- 
hole residual interaction and that g' cannot be kept equal 
to 1/3 [13, [l3|. However, for the 0+ — > 0'*" channel in 
the present paper, the direct contributions from the pion 
vanish. Therefore, in this sense, the self-consistency is 
also fulfilled in RH+RPA approach for the superallowed 
Fermi /3 decays. 



III. RESULTS AND DISCUSSION 

For all the calculations in this paper, the spherical sym- 
metry is assumed and the filling approximation is applied 



to the last partially occupied orbital. The radial Dirac 
equations are solved in coordinate space by the Runge- 
Kutta method within a spherical box with a box radius 
i? = 15 fm and a mesh size = 0.1 fm [2^. The single- 
particle wave functions thus obtained are used to con- 
struct the RPA matrices A and B in Eqs. ^ with the 
single-particle energy truncation [— A/, A/+ 120 MeV], 
i.e., the occupied states are the positive energy states 
below the Fermi surface, whereas the unoccupied states 
can be either positive energy states above the Fermi sur- 
face or bound negative energy states With these 
numerical inputs, the model-independent sum rule. 



(14) 



can be fulfilled up to 10~^ accuracy, and the isospin 
symmetry-breaking corrections Sc are stable with respect 
to these numerical inputs at the same level of accuracy. 



A. Isospin symmetry-breaking correction Sc 

In Table |T1 the isospin symmetry-breaking corrections 
6c in Eq. ([3]) for the 0+ 0+ superallowed transitions 
are shown. The results are obtained by self-consistent 
RHF-hRPA calculations with PKOl [iBl, PK02 
PK03 [i^ effective interactions, as well as by self- 
consistent RH-f-RPA calculations with DD-MEl [H], 
DD-ME2 [13, NL3 [13], TMl [H efi'ective interactions. 
The results obtained by shell model calculations (T&H) 
Q are also listed for comparison. The present corrections 
6c range from about 0.1% for the lightest nucleus to 

jabout 1.2% for the heaviest nucleus ''"'Rb, which are 2-3 
times smaller than the T&H results. It is noticed that 
even smaller values of 6c compared to the shell model cal- 
2 culations have been recently obtained in Ref. [3l[ using 

.perturbation theory. In addition, in Table Hi] the exci- 
tation energies for the 0^ 0^ superallowed tran- 
sitions corresponding to PKOl and DD-ME2 are shown 
as examples. These energies are measured by taking the 
ground state of the corresponding even-even nuclei as ref- 
erence. In the comparison with the experimental values 
taken from the recent survey [l^ , the corrections due to 
the proton-neutron mass difference in particle-hole con- 
figurations are made for the calculated results. A good 
agreement between the data and the calculated ones can 
be seen in Table HT] 

In Table U it is found that the present isospin 
symmetry-breaking corrections 6c for each nucleus can 
be unambiguously divided into two categories, those ob- 
tained by RHF+RPA calculations and those obtained 
by RH-I-RPA calculations. Comparing these two cate- 
gories, it is seen that the corrections 6c of RHF+RPA 
are systematically smaller than those of RH+RPA. On 
the other hand, it is also found that within one cate- 
gory the corrections 6c are not sensitive to specific ef- 
fective interactions or the structure of the Lagrangian 
density. For instance, within the RH+RPA framework. 
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TABLE I: Isospin symmetry-breakin g co rrections Sc for the 0^ 0^ superallowed transitions obtained by self-consistent 
RHF-I-RPA calculations with PKOl [il, PK02 [2I, and PK03 ^ as well as self-consistent RH-fRPA calculations with 
DD-MEl [111, DD-ME2 NL3 [H, and TMl gj. The column PKOl* presents the results obtained with PKOl without 
the Coulomb exchange (Fock) term. The results obtained by shell model calculations are listed in the column T&H for 
comparison. All values are expressed in percents. 





PKOl 


PK02 


PK03 


PKOl* 


DD-MEl 


DD-ME2 


NL3 


TMl 


T&H [9] 




0.082 


0.083 


0.088 


0.148 


0.149 


0.150 


0.124 


0.133 


0.175(18) 


"0 "n 


0.114 


0.134 


0.110 


0.178 


0.189 


0.197 


0.181 


0.159 


0.330(25) 




0.270 


0.277 


0.288 


0.357 


0.424 


0.430 


0.344 


0.373 


0.565(39) 


26sj _^ 26^1 


0.176 


0.176 


0.184 


0.246 


0.252 


0.252 


0.213 


0.226 


0.435(27) 


on.-i '^n-r-, 
oUo oUp 


0.497 


0.550 


0.507 


0.625 


0.612 


0.633 


0.551 


0.648 


0.855(28) 


3^Ar ^ ^'^Cl 


0.268 


0.281 


0.267 


0.359 


0.368 


0.376 


0.438 


0.320 


0.665(56) 


3»Ca 3«K 


0.313 


0.330 


0.313 


0.406 


0.431 


0.441 


0.390 


0.572 


0.765(71) 


42Ti -> -i^Sc 


0.384 


0.387 


0.390 


0.460 


0.515 


0.523 


0.436 


0.443 


0.935(78) 


26^1 _> 26^g 


0.139 


0.138 


0.144 


0.193 


0.198 


0.198 


0.172 


0.179 


0.310(18) 




0.234 


0.242 


0.231 


0.298 


0.302 


0.307 


0.289 


0.267 


0.650(46) 


3»K 3«Ar 


0.278 


0.290 


0.276 


0.344 


0.363 


0.371 


0.334 


0.484 


0.655(59) 


42Sc ^ '12^^ 


0.333 


0.334 


0.336 


0.395 


0.442 


0.448 


0.377 


0.383 


0.665(56) 


54 Co ^ 54 pg 


0.319 


0.317 


0.321 


0.392 


0.395 


0.393 


0.355 


0.368 


0.770(67) 


^•^As ^ ®^Ge 


0.475 


0.475 


0.469 


0.571 


0.568 


0.572 


0.560 


0.524 


1.56(40) 


TOgr _> ^oge 


1.140 


1.118 


1.107 


1.234 


1.232 


1.268 


1.230 


1.226 


1.60(25) 


■^"Rb ^"Kr 


1.088 


1.091 


1.071 


1.230 


1.233 


1.258 


1.191 


1.234 


1.63(31) 



TABLE II: Excitation energies for the 0"*" — > 0"*" super- 
allowed transitions measured by taking the ground state of 
the corresponding even-even nuclei as reference. In the com- 
parison with the experimental values taken from the recent 
survey [l^, the corrections due to the proton-neutron mass 
difference in particle-hole configurations are made for the cal- 
culated results. All units arc in MeV. 





expt. 


PKOl 


PKOl* 


DD-ME2 




-1.908 


-1.698 


-2.307 


-2.236 


140 ^ "N 


-2.831 


-2.420 


-2.989 


-3.081 


iBNe ^ i»F 


-3.402 


-3.195 


-3.497 


-3.451 


26si ^ 26^1 


-4.842 


-4.531 


-5.139 


-5.110 


30g _^ 30p 


-5.460 


-4.845 


-5.326 


-5.395 


34 Ar -> 34C1 


-6.063 


-5.559 


-6.129 


-6.278 


38Ca ^ 


-6.612 


-6.035 


-6.611 


-6.775 


42Ti _ 42sc 


-7.000 


-6.661 


-6.970 


-6.964 


26 Al _ 26^g 


4.233 


3.908 


4.372 


4.350 


3*01 34s 


5.492 


5.062 


5.428 


5.561 


38k ^ 38 Ar 


6.044 


5.557 


5.936 


6.083 


*2Sc -> -i^Ca 


6.426 


6.118 


6.333 


6.333 


54 Co ^ 54 pg 


8.244 


7.720 


8.221 


8.240 


66 As ^ 66cg 


9.579 


9.044 


9.488 


9.677 




9.970 


9.632 


9.805 


9.852 


^%b ^ ^"Kr 


10.417 


10.005 


10.349 


10.437 



both the Lagrangian densities with density-dependent 
mcson-nucleon couplings (DD-MEl, DD-ME2) or with 
non-hncar meson couphngs (NL3, TMl) lead to quite 
similar results. 

To understand this systematic discrepancy between 
RHF+RPA and RH-f-RPA, it must be kept in mind that 
in RHF+RPA the exchange (Fock) terms of mesons and 
photon arc kept in both the mean field and RPA lev- 



els, whereas they are neglected altogether in RH+RPA. 
Among all the Fock terms, we expect, in particular, the 
exchange terms of the Coulomb field to play an important 
role due to the following reason. The IAS would be de- 
generate with its isobaric multiplct partners, i.e., = 0, 
and it would contain 100% of the model-independent 
sum rule (|14p . i.e., Sc = 0, if the nuclear Hamiltonian 
commutes with the isospin raising and lowering opera- 
tors T±. This would be true when the Coulomb field 
is switched off. While this degeneracy is broken by the 
mean field approximation, no matter the exchange terms 
of mesons are included or not, it can be restored by the 
RPA calculations as long as the RPA calculations are 
self-consistent [13 • Therefore, the Coulomb field is es- 
sential for the O"*" 0'^ superallowed transitions and the 
Coulomb exchange (Fock) term should be responsible for 
the the different isospin symmetry-breaking corrections 
Sc in RHF-I-RPA and RH-I-RPA approaches. 

In order to verify the above argument, we have per- 
formed the following calculations. Using PKOl, the 
Hartree-Fock calculations are performed by switching off 
the exchange contributions of the Coulomb field. From 
the single-particle spectra thus obtained, self-consistent 
RPA calculations are then performed. One may notice 
that in such calculations some nuclear properties includ- 
ing binding energies and rms radii can no longer be re- 
produced. However, this docs not hinder us from dis- 
cussing the physics we are concerned with. The isospin 
symmetry-breaking corrections Sc and the excitation en- 
ergies Ex thus obtained are listed in the column denoted 
as PKOl* in Tables [T] and [TTl It is seen that these re- 
sults are almost the same as those of RH+RPA calcula- 
tions with DD-MEl, DD-ME2, NL3, and TMl, i.e., by 
switching off the exchange contributions of the Coulomb 
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field, Ex and 5c in RHF+RPA calculations recover the re- 
sults in RH+RPA calculations. In other words, although 
the meson exchange terms can be somehow effectively in- 
cluded by adjusting the parameters in the direct terms, 
this has not been done for the Coulomb part in the usual 
RH approximation. 

Therefore, one can conclude that the proper treat- 
ments of the Coulomb field is very important to extract 
the isospin symmetry-breaking corrections Sc- 



B. Nucleus-independent J^t values 

Among the 0+ 0+ superallowed transitions listed 
in Table |T1 some of their measured ft values are 
summarized in a recent survey [l9l |. To obtain the 
nucleus-independent J-t values from each experimental 
ft value, apart from the isospin symmetry-breaking cor- 
rections Sc in Table U one still needs the values of 
the transition-dependent radiative corrections and 
nuclear-structure-dependent radiative corrections Sjq^. 

Using the S'j^ and iJns values from recent calculations 
0, Sc in Table HI and measured ft values [3, the 
nucleus-independent J^t values for superallowed Fermi /3 
decays are listed in Tablc HlIl together with the average !Ft 
values and the values of chi-square per degree of freedom 
'X^ jv^ in which the uncertainty of Sc is taken as zero. 

It is found that the chi-square per degree of freedom 
X^/i^ is 1.0 '--^ 1.1 s for all effective interactions employed. 
This indicates the constancy of the nucleus-independent 
Tt values is satisfied, even though not as well as the shell 
model calculations in Ref. [l^. It is also found that the 
Tt values of RHF-I-RPA are about 2 s larger than those 
of RHH-RPA, which are larger than the difference due to 
the different effective interactions in either RHF or RH 
approximations. 

The results of RHF+RPA with PKOl, RH+RPA with 
DD-ME2 and NL3 arc plotted as a function of the charge 
Z for the daughter nucleus in Fig.[l]to illustrate the con- 
stancy of the nucleus- independent Tt values. The shaded 
horizontal band gives the standard deviation, which com- 
bines the statistical errors and jv, around the average 
Tt value. 

In order to get a deeper understanding on the treat- 
ment of the Coulomb field, the Tt values from RPA cal- 
culations using Skyrme Hartree-Fock (SHF) with SGII 
effective interaction are shown in panel (b) of Fig. [l] in 
which the isospin symmetry-breaking corrections Sc are 
taken from the Table I in Ref. [l3|. It should be empha- 
sized that in these results the exchange contributions to 
the Coulomb mean field are treated in Slater approxima- 
tion. Although this model leads to a similar average Tt 
value, Tt ~ 3081.1(7) s, it is found that the chi-square 
per degree of freedom 'x^ jv = 1.5, i.e., the constancy of 
the Tt values in this SHF framework is not as good as 
that given by the relativistic calculations. In particular, 
the Tt value deduced from the nucleus ^"'Rb is seriously 
overestimated. 



C. CKM matrix 

With the nucleus-independent Tt value, the element 
Vud of the CKM matrix can be calculated by 

"'^ 2G|,(l + A]^)Ji' ^ ^ 

where K/{hcf = 8120.2787(11) x 10"!° GeV^s, 
Gp/ihcf = 1.16637(1) X 10-5 GeV"^ 0, and = 
2.361(38)% Q. Then in combination with the other two 
CKM matrix elements = 0.2255(19) and |K&| = 

0.00393(36) 01, one can test the unitarity of the matrix. 

The element Vud as well as the sum of squared top-row 
elements of the CKM matrix are listed in Table HVl The 
uncertainties of the present results are underestimated 
to some extent as the uncertainty of Sc is assumed to be 
zero and the systematic errors are not taken into account. 
In Fig. [2l the element Vud of the CKM matrix obtained 
by RHF+RPA calculations with PKOl and by RH+RPA 
calculations with DD-ME2 are shown in comparison with 
those in the shell model (H&T) [l^ as well as in neutron 
decay Q , pion f3 decay Q , and nuclear mirror transitions 
i- 

It can be clearly seen in Table IIVI that the matrix 
element \Vud\ determined by the 0+ 0+ superal- 
lowed transitions mainly depends on the treatment of 
the Coulomb field and less sensitive to the particular ef- 
fective interactions. Switching either on or off the ex- 
change contributions of the Coulomb field, the discrep- 
ancy caused by different effective interactions is much 
smaller than the statistic deviation. It is interesting to 
note that the present values well agree with those 

obtained in neutron decay, pion (3 decay and nuclear mir- 
ror transitions. However, the sum of squared top-row 
elements considerably deviates from the unitarity con- 
dition, which is in contradiction with the conclusion in 
shell model calculations (H&T) [l^. This calls for more 
intensive investigations in the future. For example, mean 
field and RPA calculations including the proper neutron- 
proton mass difference, isoscalar and isovector pairing, 
and deformation should be done. It should also be em- 
phasized that apart from the proper treatment of pairing 
by either BCS or Bogoliubov approaches, the particle 
number projection must be implemented as well in order 
to remove the artificial isospin symmetry breaking effects 
due to the particle number violation. 



IV. SUMMARY AND PERSPECTIVES 

In summary, self-consistent relativistic RPA ap- 
proaches are applied to calculate the isospin symmetry- 
breaking corrections Sc for the 0+ 0+ superal- 
lowed transitions. In the RHF+RPA framework the 
density-dependent effective interactions PKOl, PK02, 
and PK03 are employed, while in the RH+RPA frame- 
work the density-dependent effective interactions DD- 
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TABLE III: Nucleus-independent Tt values. The average Tt value and the normalized /v appear at the bottom. All units 
are in s. 







PKOl 


PK02 


PK03 


PKOl* 


DD-MEl 


DD-ME2 


NL3 


TMl 






3079.6(45) 


3079.5(45) 


3079.4(45) 


3077.5(45) 


3077.5(45) 


3077.5(45) 


3078.3(45) 


3078.0(45) 


"0 -> 


"N 


3078.2(31) 


3077.5(31) 


3078.3(31) 


3076.2(31) 


3075.8(31) 


3075.6(31) 


3076.1(31) 


3076.8(31) 


34 Ar^ 


34 CI 


3081.9(84) 


3081.5(84) 


3082.0(84) 


3079.1(84) 


3078.8(84) 


3078.6(84) 


3076.7(83) 


3080.3(84) 


26 Al^ 


2'^Mg 


3077.7(13) 


3077.7(13) 


3077.5(13) 


3076.0(13) 


3075.8(13) 


3075.8(13) 


3076.6(13) 


3076.4(13) 


34Cl-> 


34g 


3083.5(16) 


3083.3(16) 


3083.6(16) 


3081.6(16) 


3081.4(16) 


3081.3(16) 


3081.8(16) 


3082.5(16) 




38 Ar 


3084.1(16) 


3083.8(16) 


3084.2(16) 


3082.1(16) 


3081.5(16) 


3081.3(16) 


3082.4(16) 


3077.8(16) 




42Ca 


3082.7(21) 


3082.6(21) 


3082.6(21) 


3080.7(21) 


3079.3(21) 


3079.1(21) 


3081.3(21) 


3081.1(21) 


54Co^ 


^4Fe 


3083.9(27) 


3083.9(27) 


3083.8(27) 


3081.6(27) 


3081.5(27) 


3081.6(27) 


3082.7(27) 


3082.4(27) 




^4Kr 


3094.8(87) 


3094.7(87) 


3095.3(87) 


3090.3(87) 


3090.2(87) 


3089.4(87) 


3091.5(87) 


3090.2(87) 


average 




3081.4(7) 


3081.3(7) 


3081.4(7) 


3079.5(7) 


3079.1(7) 


3079.0(7) 


3080.0(7) 


3079.1(7) 






1.1 


1.1 


1.1 


1.0 


1.0 


1.0 


1.0 


1.0 





3120 


" PKOl 








SGII 


' 1 ' 


1 i i 1 








3100 












































3080 














i I 










(a) : 






(b) : 
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FIG. 1: Nucleus-independent values as a function of the charge Z for the daughter nucleus. The values of 5c are respectively 
obtained by RHF-fRPA calculations with PKOl (a), by RH-I-RPA calculations with DD-ME2 (c) and NL3 (d), as well as by 
SHF+RPA calculations with SGII [l^l (b). The shaded horizontal band gives one standard deviation around the average n 
value. 



TABLE IV: The element Vud and the sum of squared top-row 
elements of the CKM matrix. 







|Kd|' + |K,.|' + |K,6|' 


PKOl 


0.97273(27) 


0.9971(10) 


PK02 


0.97275(27) 


0.9971(10) 


PK03 


0.97273(27) 


0.9971(10) 


PKOl* 


0.97303(26) 


0.9977(10) 


DD-MEl 


0.97309(26) 


0.9978(10) 


DD-ME2 


0.97311(26) 


0.9978(10) 


NL3 


0.97295(26) 


0.9975(10) 


TMl 


0.97309(26) 


0.9978(10) 



MEl and DD-ME2 as well as the nonlinear effective in- 
teractions NL3 and TMl are used. 

It is found that the proper treatments of the Coulomb 
field is very important to extract the isospin symmetry- 



breaking corrections Sc- By switching off the ex- 
change contributions of the Coulomb field, and 6c in 
RHF-fRPA calculations recover the results in RH-f RPA 
calculations. In other words, although the meson ex- 
change terms can be somehow effectively included by ad- 
justing the parameters in the direct terms, this has not 
been done for the Coulomb part in the usual RH approx- 
imation. 

With the isospin symmetry-breaking corrections 6c 
calculated by relativistic RPA approaches, the nucleus- 
independent J^t values are obtained in combination with 
the experimental ft values in the most recent survey 
and the improved radiative corrections. It is found that 
the constancy of the !Ft values is satisfied for all self- 
consistent relativistic RPA calculations here. It is also 
found that the values of RHF+RPA are about 2 s 
larger than those of RH+RPA, which arc larger than the 
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nuclear 

neutron 

superallowed ^^^^^ pion p 
transitions decay 



/ I \ 



i 

H&T 



PK01 



DD-ME2 



nuclear 
mirror 
transitions 



ther RHF or RH approximations. 

The values of \ Vud\ thus obtained weU agree with those 
obtained in neutron decay, pion (3 decay, and nuclear 
mirror transitions. However, the sum of squared top-row 
elements considerably deviates from the unitarity con- 
dition, which is in contradiction with the conclusion in 
shell model calculations (H&T) p^ . 

For the further studies, more intensive investigations 
including the proper neutron-proton mass difference, 
isoscalar and isovector pairing, and deformation should 
be done. 



FIG. 2: The element Vud of the CKM matrix obtained by 
RHF+RPA calculations with PKOl and by RH-I-RPA calcu- 
lations with DD-ME2 in comparison with those in shell model 
(H&T) [l§] as well as in neutron decay 01 , pion /3 decay 0) 
and nuclear mirror transitions [a|. 



difference due to the different effective interactions in ci- 
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